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Abstract 

Graphene outstanding properties directly come from its pecular electronic structure and thus 
from the honeycomb lattice symmetry. The way interaction with the substrate impact this lattice is 
of primary importance. This is peculiarly true for epitaxial graphene because of the SiC substrate. 
The advantage of this system that produces macroscopic samples in registry with a substrate could 
turn to a major drawback if the graphene lattice reveals to be strongly distorted. Extensive ab initio 
calculations supported by Scanning Tunneling Microscopy experiments, demonstrate here that the 
substrate indeed induces a strong nanostructuration of the interface carbon layer. It propagates 
to the above C layer where it generates incommensurate ripples in the honeycomb lattice. 
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I. INTRODUCTION 



Between 2004 and 2005, three groups propelled graphene to the foreground 
a. They brought .t from the theorefcal to the real experimental world by pfrrthetizing it 

for the first time. Indeed, since the pioneering work of Wallace in 1947 [5] very exotic 
properties had been forcast for this material that remained untested. 3 years after its first 
synthesis, graphene is up to its promises : room temperature quantum Hall effect has been 
evidenced [6[ together with large coherence length and high electronic mobility 
Furthermore, these fundamental properties give a high application potential to graphene, 
especially in nanoelectronics (qbits, transistors) [Sj. Outstanding electronic properties 
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ll| together with a high potential for application explain why graphene is becoming the 



new star of condensed matter physics. 

Graphene is the name given to one isolated plane of carbon atoms arranged on a honey- 
comb lattice. The unit cell contains two equivalent carbon atoms (named A and B). The 
interaction with this pecular lattice makes electrons behave like relativistic massless fermions 
governed by a Dirac like equation. Indeed the band structure of graphene shows two bands 
with a linear dispersion that cross at the Fermi level. The crossing point is called the Dirac 
point. The anomalous quantum Hall effect of graphene for example is a direct consequence 
of this band structure. 

Two ways are used to synthetize high crystalline quality samples : either from mechanical 



exfoliation of highly oriented p yro litic graphite (HOPG) 



from 4H or 6H SiC surface 



!)vrontic g 
MM- 



or from sublimation of Si 



We focus here on epitaxial graphene on SiC since 



then graphene is directly grown on a substrate and this approach is the only one that could 



been s 
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lown to produce from few 
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16|. 



lead to mass production for nanoelectronics needs. It has 
layers of graphite (FLG) to one unique graphene layer 0, 

High quality graphene is required for transport properties and the development of nano- 
electronic applications . The mastering of the interface morphology is then a priority. The 
main question one has to answer is how the SiC substrate impacts the grap hene unique struc- 



ture. Indeed, Angle Resolved Photoelectron Spectroscopy (ARPES) [16|, [iT], ll8|, 13, 2ol. |21| 
and transport measurements [3] showed that the system morphology has a strong effect 
on the electronic structure of the FLG. The crystalline structure of the system has been 
approached by different techniques (x-ray diffraction 22l, l23|, l2J], Low Energy Electron 
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neling microscopy (STM) 
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Diffraction (LEED) [251, l2aj27Lcore level shift (CIS) (l8|, l20|, l25|, l26| and scanning tun- 
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28j). The existence of a QV3xQV3R30 (hereafter 6R3) 



cell common to the SiC substrate and the graphene is rather well established for the Si- 
terminated surface. The C-terminated case is not so clear but the interface geometry does 
not seem to involve the 6R3 cell and calculations with this geometry are given here only for 
direct comparison. 

Ab initio calculations using a highly simplified geometry for the interface (a y^x v^-R30 
cell, hereafter R3) have established the bufferlayer role of the first C layer that electronically 



30[. Here we address the actual 



decouples the graphene planes from the substrate 
6v^x6v^-R30 cell and on the basis of extensive ab initio calculations supported by scanning 
tunneling experiments we demonstrate that the substrate has a strong effect on the first two 
carbon layers. The interface C buffer layer is strongly distorted so that it cannot be seen 
as a graphene layer but it should instead be considered as the very beginnings of graphene. 
Its actual structure is elucidated at the atomic scale. Moreover, ripples are shown to be 
generated in the otherwise graphene like second plane. 

Calculation and experimental details are given in part II. The third part presents the ab 
initio results for both Si- and C-terminated faces, for the buffer layer and for the first actual 



graphene layer. STM results are compared to ab initio calculations in part IV. 



II. CALCULATION AND EXPERIMENTAL DETAILS 



Ab initio calculations have been performed whithin the Density Functionnal Theory using 
the code VASP [31]. Ultra soft pseudopotentials (USPP) |32| have been used with a plane 
wave basis cutoff equal to 211 eV; The USPP have been extensively tested, especially, the C- 
short USPP was shown to correctly reproduce the band structure of graphene, graphite and 
bulk and surfaces of different SiC polytypes 29 1. The Perdew and Wang [sll formulation of 
the General Gradient Approximation is used. Brillouin zone integration is performed using 
the r point. 

Calculations have been performed with one or two carbon layers on the Si-terminated 
and C-terminated surfaces of SiC. The graphene and the (0001) SiC lattices are nearly 
commensurate with a 6R3 common cell with respect to SiC or a 13x13 cell with respect 
to graphene. In the following, periodicity expressed with respect to SiC (graphene) will 
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be referred to -SiC (-G) for clarity. The calculations are performed in the actual 6R3- 
SiC interface geometry. We used complete and flat graphene layers as a starting point for 
all the C planes above the SiC substrate. This is to account for the ARPES results that 
demonstrated the existence of the a bands skeleton for the interface C layer j2(]]: these 
bands are related to the presence of a well developped honeycomb like lattice. When an 
additional graphene layer is added, the converged geometry of the interface C layer remains 
unchanged and forms the so-called bufferlayer. 

The starting buffer layer system (1 C layer on SiC) was described by a supercell formed 
of 4 SiC bilayers plus one layer of C atoms arranged on a honeycomb lattice. H atoms are 
bound to the second SiC surface to passivate it and the whole cell contains 1310 atoms. 
The atomic positions where relaxed so that residual forces on the atoms are lower than 
0.015 eV/A. The same cell with an additional graphene layer would be too large to keep the 
computational time within a reasonnable value (several months had already been necessary 
to reach convergence for each system) and we had to restrict the cell to two SiC bilayers, 
the buffer layer and the graphene layer (1216 atoms). Atoms in the lower SiC bilayer and 
the subsequent C plane are kept fixed to bulk positions. Because of the limited number of 
planes used to describe the substrate, forces could not be zeroed in these 3 layers that remain 
constrained. However, forces in the above graphene, buffer layer and last SiC plane could 
be relaxed to negligible values and the calculation that started from two fiat honeycomb 
lattices reproduces the buffer layer geometry found with the first cell. 

STM experiments were performed at room temperature, in ultra-high vacuum, using 
a home made microscope. Graphitization of n-type (nitrogen 1x10^^ cm~^) 6H-SiC(0001) 
substrates was achieved by successive annealing into UHV, controlled by LEED and Auger 
spectroscopy. Experiments were done first on the buffer layer surface, and second on a mixed 
surface with mono- and bilayer graphene areas. The details of the growth process and the 
determination of layer thickness have been discussed elsewhere [l^. We used the symmetric 
STM atomic contrast to identify areas corresponding to monolayer graphene. Results have 
been checked on 2 different substrates with different tips and at different temperatures. 
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III. AB INITIO RESULTS 



A. Buffer layer on the Si-terminated SiC surface 

Fig. la shows a large scale image of the ab initio total charge density of an uncovered 
buffer layer (the picture is identical when it is covered by a graphene layer). It has three 
main characteristics: i) an obvious apparent 6x6-SiC modulation ii) The 6R3-SiC periodicity 
iii)low regions, separated by boundaries, that form nanograins with an unexpected local 2x2- 
G symmetry. This gives a mosaic structure to the buffer layer. 

The 6x6-SiC modulation is strongly apparent (cell in red in Fig. la) because it is related 
to the bright spots of Fig. la and Fig. lb. It is not a real reconstruction since it is not 
compatible with translationnal symmetry within the buffer layer: Fig. la and Fig. lb show 
that the local environment of the atoms at the 6x6-SiC bright spots is not the same. The 
6R3-SiC periodicity (cell in blue in Fig. la) corresponds to the SiC and graphene common 
cell. It is imposed by the calculation. 

The mosaic pattern is composed of grains of more or less irregular hexagonal shape, 
20Awide. The origin of the mosaic structure comes from the superposition of the C hon- 
eycomb and SiC lattices [25^]. The two lattices do not adjust to each other. The grains 
(dark area) corresponds to regions where the SiC and honeycomb lattices match (in a local 
2x2-G or R3-SiC symmetry). Si-C bonds are formed there. It generates the small hexagonal 
patterns of Fig. la. The local 2x2-G patterns are shifted in adjacent grains, the boundaries 
also accomodate for this shift. Atoms that are not in registry lay higher above the substrate 
and form boundaries (light area). 

The Fourier transform of the charge density map of the free buffer layer (Fig. Id) and of 
the buffer layer covered with a graphene layer (not shown) are identical. Apart from the 
Ixl-G graphene and the 6x6-SiC spots, one remarkable feature of this FT image is the rather 
intense spots located around midway from the Ixl-G spots (but off-axis). These spots are 
also found in the FT of the STM images as shown below. They belong to the reciprocal 
lattice of the 6R3-SiC. Their intensity is locally enhanced (in k-space) by the presence of 
small grains with local 2x2-G (or R3 -SiC) symmetry. We call them 2x2-G spots in the 
following. Additionally, one remarks that the first order spots of the 6R3 reconstruction are 



missing at the centre of Fig. Id (although higher order spots of the 6R3 are present) [34 1. 
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The corrugation of the buffer layer is rather large with a lower to upper atom height 
difference close to 0.12 nm (Fig. If). In contrast to previous models proposed for the buffer 



layer 25|, |28| , the present structure directly comes from atomic relaxation during ab initio 



calculation starting from a flat carbon honeycomb lattice on top of the SiC surface in the 
actual 6R3-SiC common cell. The resulting morphology is in agreement with CLS data that 
evidence the existence of 2 types of C atoms that are bound or not to Si atoms 
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Ab initio calculations in the 2x2-G structure [29|] demonstrated that the interface C layer 
has not the electronic structure of graphene but acts as a bufferlayer since it allows growth of 
subsequent C planes with graphene like dispersion . The present calculations, performed in 
the actual interface geometry, support this results. It demonstrates that part of the interface 
layer indeed has the simple 2x2-G geometry. It also further evidences the buffer role of the 
first C layer since a second C layer, on top of it, presents a (wavy) honeycomb lattice. 



B. Graphene layer on the Si-terminated SiC surface 

The nanostructuration induced by the substrate in the bufferlayer propagates to the ontop 
graphene layer. It generates ripples that show up as a 6x6 modulation of the honeycomb 
lattice in the ab initio total charge density (Fig. 2a, b and c). This is confirmed by the 
Fourier transform that exhibits spots related to Ixl-G and 6x6-SiC periodicity (Fig.2d). 
High regions of the buffer layer correspond to high regions of the graphene layer (Fig. 2c). The 
ripples have an amplitude of 0.04 nm for a wavelength of 1.9 nm. Such a long wavelength 
modulation does not discriminate between the two basis atoms, A and B, of the honeycomb 
lattice. So that no A versus B contrast can be seen on the ab initio total charge density 
map in agreement with STM image (inset Fig.4b). 

An ab initio calculation of a free rippled graphene layer, where the atomic positions 
are frozen in the presently calculated 6x6 positions, reveals an ideal graphene like linear 
dispersion. The induced distorsion differenciates areas of the graphene layer but it does 
not differenciate on the whole A atoms from B atoms. The band structure calculation of 
the graphene layer with the substrate and buffer layer is unfortunately out of reach for the 
moment (more SiC bilayers need to be considered and the system would be too large for 
reasonable computational time). However, if a gap is opened in the graphene band structure 
because of the substrate, the phenomenon involved is more complex than a simple break 
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of AB sublattice symmetry as it has been proposed recently 2l|. Furthermore, ARPES 
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20| , STM IJ, |25|] and prehminary calculations 
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30j demonstrate that the buffer 



layer has no vr states in the vicinity of the Fermi level. So these states cannot be involved in 
the disymmetrization of the A and B sublattices even if the buffer and the graphene layer 
stacking is graphite like (Bernal stacking). 



C. The C-terminated case 



The total charge density maps for the C-terminated interface are given for compari- 
son with the Si-terminated interface. They give indications on what could happen at the 
C-terminated interface if the interface geometries were identical. In fact, the actual C- 
terminated geometry is not precisely known. The supercell used for these calculations is 
equivalent to the one used for the Si-terminated interface. The same complex mosaic struc- 
ture appears for the first carbon layer (Fig. 3a) but the hexagonal shape of the grains seems 
more regular and more pronounced. All boundaries intersections give rise to bright spots 
so that the 6x6 modulation is not so apparent. By the way, the nanostructuration of the 
graphene layer also appears to be stronger (Fig.Sb). The difference between the two faces 
might come from the strength of the substrate - C layer bonds. In the C-terminated case, 
covalent C-C bonds are formed. They are stronger and shorter than Si-C bonds. The in- 
duced buffer layer lattice distorsions and subsequently, the graphene layer structuration are 
stronger. 



IV. STM RESULTS ON THE SI-TERMINATED SURFACE 

The pecularities of the interface evidenced by ab initio calculations are also found in 
Constant Current topographic STM images of the buffer layer (Fig. 4a), of the first graphene 
monolayer (Fig.4b) and of their respective FT (Fig.4c and 4d). The fact that STM images are 
related to the surface local density of states whereas total charge densities are obtained from 
ab initio calculations only has a minor effect in the present discussion. For the uncovered 
buffer layer, (the so-called nanomesh or 6R3 phase in the litterature [251, |27[), a honeycomb 
pattern with the 6x6-SiC periodicity is experimentaly observed (Fig.4a), with hexagons of 
irregular shapes in accordance with the calculation (Fig. la). Although atomic resolution 
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is hardly achieved on this surface Q, 27 1, the FT agrees with the 2x2-G periodicity 
calculated for the grains. Indeed, we find some tiny spots (circled in Fig. 4c) at positions 
close to half the reciprocal vector of the Ixl-G lattice (k=29.5 nm^^). The same spots are 
found on the FT of the corresponding calculated total charge shown in Fig. Id. 

STM images of the first graphene layer capping the buffer layer strongly support our 
calculated results. A survey of the image of Fig. 4b, recorded at sample bias -0.2V, reveals 
both the graphene Ixl-G honeycomb pattern superimposed to a 6x6-SiC superstructure. 
These features are recovered in the FT shown in Fig. 4d. A careful study of STM images 
of a few nm wide, like the 2.5 x 2.5 nrn? image shown in inset, agrees with previous results 
of [l^: we do not find any detectable asymmetry between the adjacent carbon atoms of 
the graphene layer, and we use the symmetric contrast as a signature of the monolayer (an 
asymmetric contrast is found for graphene bilayer in AB stacking). As discussed in section 
JIB, we believe that this is in contradiction with the mechanism suggested by Zhou et al 



2l| for a possible gap opening at the Dirac point. 



We can take advantage of the STM sensitivity to the interface states [24, l35|, l36[ to get 
new insights into the buffer layer structure capped by the graphene top-layer. At low sample 
bias (here -0.2V), the tip probes the vr-like states of graphene, together with electronic states 
lying within the buffer layer or at the interface with SiC p^]. In the present case, interface 
states give rise to additional features in the FT of the STM image related to the buffer 
layer. As shown in Fig. 3d, we find a clear signature of the local 2x2-G symmetry, arising 
as pronounced bright spots (circle) close to the 2x2-G expected value. They correspond to 
the spots circled in the FT of the total charge density of the free buffer layer (Fig. Id). Both 
theory and STM experiment show that the buffer layer geometry is essentially unchanged 
when it is buried below one graphene monolayer. 

From this results, it appears that the C buffer layer is not an intermediate sublimation 
phase but remains at the interface. A possible scenario for growth might proceed as follows: 
Si atoms sublimates, the C atoms rearrange and form the buffer layer. More Si atoms 
disappear, the C atoms freed form a new buffer layer while the former one flattens to 
form a graphene plane. Since the buffer layer needs to form prior to a new graphene layer 



appearance, a layer by layer growth is favoured as observed experimentally ly, |24| |. 
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V. CONCLUSION 



As a conclusion, the C buffer layer is not an intermediate sublimation phase. It is always 
present at the Si-terminated SiC - graphene interface and it decouples the graphene layer 
from the substrate. Most of the system interest relies on the presence of this layer. It has 
a mosaic structure that derives from a honeycomb lattice distorted because of C-Si bond 
formation. Its structure is complex, it shows an 6x6-SiC modulation also found in STM 
images while the actual common cell corresponds to a 6\/3x6y/3R3Q-SiC. It also involves a 
local 2x2 symmetry (with respect to graphene) in the region where the honeycomb and the 
SiC lattices nearly match. This unexpected local symmetry appears in the Fourier transforms 
of the STM images confirming the existence of matching zones. The existence of this buffer 
layer agrees with a layer by layer growth mechanism where a new buffer layer is formed 
at the interface while the former one evolves into a graphene layer. The substrate induced 
deformation propagates to the first graphene layer where it creates incommensurate 6x6-SiC 
ripples. This incommensurate modulation, observed on STM images too does not break the 
AB symmetry. The effect of the ripples on transport properties requires further investigation. 
The possible use of the graphene and buffer layer substrate induced nanostructuration to 
form arrays of nanometric size islands also deserves further consideration. 
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Figure captions 

Figure 1: (color on line) Buffer layer ontop of the Si-terminated SiC surface. The SiC 
- graphene common cell is shown in blue and the incommensurate 6x6-SiC modulation in 
red. a) Mosaic structure of the C buffer layer evidenced in this 11x11 nw? image of the 
total ab initio charge density, 0.04 nm above the upper atom, b) total charge density in 
the 6R3-SiC unit cell, c) cross section of the total charge density along the hue define in 
b) . The black dots that appear when the cross section goes through the middle of an atom 
are due to the use of pseudopotentials (no core electrons), d) Fourier transform of the large 
scale image of the total charge density a). The solid line arrow points to a Ixl-G spot, the 
dashed line arrow to a 6x6-SiC spot and the spots in the circle corresponds to the local 
2x2-G periodicity, e) positions of the atoms in the unit cell. Si atoms in the last SiC plane 
are in blue, the colour of the atoms in the buffer layer varies as a fonction of their height, 
ranging from green (close the substrate) to red (uppermost atom), f) Height profile of the 
buffer layer atoms along the line defined in e. The line is just a guide for the eyes. 

Figure 2: (color online) Graphene layer above the C buffer layer for the Si-terminated 
surface, a) large scale image (11x11 nw?) of the total charge density 0.04 nm above the 
uppermost atom in logarithmic scale. It exhibits the characteristic 6x6-SiC periodicity (cell 
shown in red). The 6R3-SiC common cell is shown in blue, b) total charge density in the 
6R3-SiC unit cell, c) cross section of the total charge density along the hue define in b). 
The section is not taken at the same position than in Fig.lc to focus on the graphene layer, 
d) Fourier transform of the large scale image of the total charge density a). The solid line 
arrow points to a Ixl-G spot, the dashed line arrow to a 6x6-SiC spot, e) positions of the 
atoms in the unit cell, C atoms in the buffer layer are in blue, the colour of the atoms in 
the graphene layer varies as a fonction of their height, ranging from green (close the buffer 
layer) to red (uppermost atom), f) Height profile of the graphene layer atoms along the hne 
define in e. 

Figure 3: (color online) Morphology of the buffer layer and the ontop graphene layer for 
the C-terminated SiC surface.a) charge density map of the interface buffer layer (8x8 nm^), 
b) charge density map of the first ontop graphene layer (8x8 nrn^). 

Figure 4: (color onhne) STM images and their related Fourier transforms of the buffer 
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layer and the graphene layer on the Si-terminated surface, a) 12x12 nm^ large STM image 
of the buffer layer (V=-2.0 V and 1=0.5 nA). It clearly shows the irregular hexagonal pattern 
of the 6x6 -SiC periodicity, b) 12x12 nm^ STM image of the graphene monolayer (V=-0.2 
V, 1=0.1 nA) (inset 2.5x2.5 nrn?). c) Fourier transform of image a), the dashed line arrow 
points to a 6x6-SiC spot and the 2x2-G spots are circled, d) Fourier transform of image 
b), the arrow points to the Ixl-G spot, the dashed arrow points to the 6x6-SiC spots and 
2x2- G related spots are circled. 
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